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Abstract 

This work studies the stability of solitary waves of a class of sixth-order Boussinesq equations. 

1 Introduction 

In this work we study the generalized sixth-order Boussinesq (GSBQ) equation [5, 8, 9] 



^" ^tt — '^xx I pUxxxx I ^xxxxxx \J\'^))xx V J 

pg where / G C^ is homogeneous of degree p > 2. Neglecting the sixth-order term, equation (1.1) becomes 

^^ a generalization of the classical Boussinesq equations 

_• utt = u^x + Purcxx^- {f{u))x^, f3 = ±l, (1.2) 

^^ Equation (1.2) was originally derived by Boussinesq [4] in his study of nonlinear, dispersive wave 

^~~^ propagation. We should remark that it was the first equation proposed in the literature to describe 

". . this kind of physical phenomena. Equation (1.2) was also used by Zakharov [24] as a model of nonlinear 

^ string and by Falk et al [11] in their study of shape-memory alloys. 

When /3 = 1, equation(1.2) is called "bad" Boussinesq equation, while (1.2) with /3 = —1, 



Utt = Uxx - Uxxxx - {f{u))xx, (1-3) 

is called "good" Boussinesq equation. Given certain conditions on /, (1.3) possesses special traveling- 
wave solutions with finite energy. Indeed, (1.3) can be written as the system of equations 

(1.4) 
ut = (u - u^x - f{u))x 

By a solitary wave solution of (1.4), we mean a traveling- wave solution of the form 0{x — ct), vanishing 
at infinity, where c is the speed of wave propagation. It was shown in [3, 17] that these solutions are 
of the form if — [ip, —cip) so that they must satisfy 

(1- c^)^ - ^" - /(^) = 0. (1.5) 
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Bona and Sachs in [3] proved that the solitary waves of (1.4) are stable under an appropriate 
convexity condition. Liu [17, 18] showed the nonlinear instability of solitary waves of (1.4). His proof 
was based on a modification of the general argument of [13]. 

Equation (1.1) can be also written as the following system of equations 

Vt = {u + (iu^^ + U^xxx - I{u))x 

If we put the solitary wave form if{x — ct) into (1.1), we obtain 

(l-c2)^ + ^¥." + ^""-/(^) = 0. (1.7) 

It is worth noting that the solitary wave solutions of equation (1.7) have been investigated numerically 
and the two classes of subsonic solutions corresponding to the sign of /3 have been obtained, more 
precisely, the monotone shapes and the shapes with oscillatory tails [5]. 
The system (1.6) has the conserved quantities 

E(u,v)^ f hul,-f3ul+u' + v')-F{u)dx (1.8) 



2 
Q{u,v) — I uv dx (1-9) 



We also note that, at least formally, the quantity 

udl'^v dx 



is conserved for any positive integer k. If i^ is a solution of the solitary wave equation (1.7), then 
(p = {ip,—ap) satisfies 

E'{0) + cQ'{0)^Q, 

so solitary waves are critical points of the action 

L{u,v)=E{u,v)+cQ{u,v). (1.10) 

Our aim here is to study the stability of solitary waves of (1.1). 

This paper is organized as follows. In Section 2, we consider the properties of ground state solitary 
wave solutions. The solitary wave equation (1.7) is a fourth-order elliptic equation, and is identical, 
after a rearrangement of parameters, to the solitary wave equation that arises in the study of the fifth- 
order KdV equation. The variational, regularity, and decay properties of this equation were considered 
in [15], so we refer to this work for several results. In Section 3 we prove the main stability result. 
Theorem 3.2, which states that the set of ground state solitary waves is stable if d"{c) > 0, where d is 
defined by equation (3.6). In Section 4 we prove the main instability result, Theorem 4.2, which states 
that a given ground state is orbitally unstable if there exists an "unstable direction". In Theorem 
4.3 we show that such an unstable direction exists provided d"{c) < 0. Using a different choice of 
unstable direction, we also derive in Theorem 4.4 explicit conditions on p, /3 and c that imply orbital 
instability. Section 5 is devoted to establishing further properties of the function d. We first show 
that when f{u) = \u\p~^u for p < 5, there exist c near c* such that d"{c) > 0. See Theorem 5.1 and 
Corollary 5.1. We then derive in Theorem 5.2 the main scaling identity satisfied by d, and use it to 
prove that d"{c) may change sign at most once along each semi-ellipse in the (/3, c)-plane. Finally, in 
Section 6, we outline the numerical method used to compute the function d, and present the results 
of these numerical calculations. The main conclusions that can be drawn from these results are found 
in Observation 6.1. 



Notations 

For each r € M, we define the translation operator by r^w = u{- + r). 

Given a sohtary wave (p of (1.6), the orbit of is defined by the set O^ — {Tr^p; r G M}. 

We shall denote by g the Fourier transform of g, defined as 

?(C) = / g{u:)e-'-< dc.. 



For s e M and 1 < p < oo, we denote by H^-p{M.), tlic Bcsscl potential space defined by H'''P{ 
A^*iP(M), with respect to the norm 



\\9\\H^.p(m = IIA^'^I 



LP( 



where A^ = (/ — 9^)*/^. In particular, we define the nonhomogeneous Sobolev space H''{R) = if'*'^(M). 
Let ^ be the space defined by 

jT = ij2(M) X l2(M), 

with the norm 

||u||^. = ||(u,w)||^ = ||u||h2(r) + ||i'|Il2(r). 

For any positive numbers a and b, the notation a < b means that there exists a positive (harmless) 
constant ^ such that a < lib- We also use a ^ b when a < b and b < a. 

2 Existence of Solitary Waves 

Solutions of the solitary wave equation (1.7) may be shown to exist via the following variational 
problem. Define 

I{u} ^ f ul^~ Pul + (1 - ^)u^ dx (2.1) 

K{u) = (p + 1) / F{u) Ax (2.2) 



where F' — f and F(0) — 0. When c^ < 1 and /? < /?* = 2-\/l — c? (equivalently when /3 < 2 and 



c^ < cl, where c* — Jl — /3^/4 and /?+ ~ max{/3,0}), the functional / is coercive in the sense that 

/(«)>C(/3,c)h||2,,(jj) (2.3) 

where 

^('^'^)>{i-c^-i;yr^ ?>o}>«- 

Since K{u) < C\\u\\'^2ib^), it follows that for A > we have 

Mx = inf{/(u) \ueH^{R), K{u) = A} > 0. 

We say that a sequence Uk is a minimizing sequence if K{uk) — > A > and I{uk) — >■ M^. The following 
result is a consequence of the concentration-compactness theorem, and was shown in [15] for a more 
general class of homogeneous nonlinearities (see also [10, 14]). 

Theorem 2.1 Fixp > 1. Suppose c^ < 1 and l3 < f3^. If Uk is a minimizing sequence for some A > 0, 
then there exists a subsequence Uk- , scalars yj and Tp & if^(M) such that Uk{- ~ i/j) ~> ip i''^ H'^{M). 



Since the function -0 achieves the minimum Mx it satisfies the Euler-Lagrange equation 

for some multipher /i. Muhiplying this equation by f/; and integrating over M, it follows that M\ = 
I{ip) = /i.(p + 1)A, so fj, > 0. Thus if = n^^^P~^^^ is a solution of the solitary wave equation (1.7). Such 
solutions are referred to as ground states and, by the homogeneity of F, achieve the minimum 



"^(/^-^) = ^"f { jr(Jffp+i) ■ " ^ ^'W'" ^ o} 



The set of all ground states will be denoted by ^(/3,c). Multiplying the solitary wave equation (1.7) 
by (fi and integrating gives I{ip) — K{ip). Thus the set of ground states is given by 

^(/3, c) = {^ e H^R) : /(^) = K{^) = m(/3, c)^}. (2.4) 

We shall denote 

^1/?, c) = {^ = (^, ~ap) e ^ : ^ e ^(/?, c)}. 

As mentioned in the introduction, elements of '^{P, c) are critical points of the action L defined by 
(1.10). In fact, elements of ^(/3,c) are minimizers of L subject to the constraint P = 0, where 

P{w)^ {L'{w),w). (2.5) 

Theorem 2.2 Suppose j3 < j3^, and c^ < 1. Let 

yK :^{we -ST; w^O, P{w)^0}. (2.6) 

The following are equivalent. 

(a) (/? e c/K and L{0) = inf{L(uJ) : w G ^/K}. 
Proof. The identities that we shall need relating the two variational problems are 

L(w, v) = \l{u) - -^K{u) + \ [ [cu + vf Ax (2.7) 

and 

P{u, v) ^ I{u) - K{u) + / (cm + vf dx. (2.8) 



From this it follows that, for any (u,w) G ^, we have L{u,v) = ^? ,^s K(u). 

First suppose G ^(/3,c). Then by definition /((/j) = K{ip), so P(i^) = and thus G ^. 
Denote A = K{(p). Then /(i^a) minimizes /(u) over all u G -ff^(M) such that K{u) = A. Now let 
w — {u,v) G ^. Then K{u) > 0, so if we set u = an where a — {K{if)/K{u))p+T^ , then K(u) — K{ip) 
and consequently I{ip) < I{u). Therefore 

= P{ip) = I{ip) - K{lp) < I{u) - K{u) = a^I(u) - aP+^K{u) ^a^{l~ aP-^)I{u), 

which implies a < 1. Thus K{ip) < K{u), and it follows that 

^^' 2{p+l) ^^' - 2{p+l) ^ ' "■ ' 



Hence (i) implies (ii). 

Next suppose = ((^, tp) £ Af solves the minimization problem. We need to show that ip G ^(/3, c) 
and ip — —c(p. Denote A — K{tp) > and suppose u e iJ^(M) minimizes / subject to the constraint 
K{-) = A. Then 

Uxxxx + fiuxx + (1 - c^)w = Ai/(u) 
for some /i. Multiplying by u and integrating gives I{u) — ^K{u) ~ fiX. Since 

I{u) < I{ip) = K{ip) - I {cip- ipf dx < K{ip) = A, (2.9) 



we have /i < 1. On the other hand, if we set ii — fip-^ u, then I{u) — K{u) so if we define w = (u, —cu) 
then we have w G J\^ . Therefore L[ip) < L{w). Since £ -jV vfe have L{0) = ^?~'J^^s K{ip) and thus 



2b +1) 

< K{w) 



2(P+1) ^ ^ 

P — 1 2 

^ 2(pTT)^^^(^) 



It then follows that /i > 1 and thus /i = 1. This implies I(u) ~ K{u) = A. But (2.9) then implies that 
I{(p) = K{ip) = A and ^ = —c(p, so we have f £ ^{(3, c) and therefore (p £ ^(/3, c). This completes the 
proof. n 

As shown in [15], solitary waves have the following regularity and decay properties. 

Theorem 2.3 Suppose ip £ iJ^(M) is a weak solution of (1.7) and that / £ C''{R). Then <p is a 
classical solution and p £ C'^''"*(M). Furthermore, p decays exponentially as \x\ — )■ oo. 

It is noteworthy that regularity and decay properties of the solutions of (1.7) can be obtained by 
using an argument similar to [10] via the following equivalent form of (1.7) 

ip = k* f{ip), 

where 

m - g._^g.Vi^,. ^ (2-10) 






Figure 1: The kernel k, shown here for c = 1/2, and /3 = —2, f3 — and /3 — 1.5. 



c2 < 1 and ^ < /3* = 2VT 
expressions for k. 



Using the residue theorem, one obtains the following explicit 



( vr f J_ 



g-Ai|a:| ]_Q->'2\x\ 



K{x) = <^ 



= < iiVl 



-V^i^i 






(2.11) 



^ 2^!,(J+c^2) (a;cos(a;a;) + crsin(^|x|)) , ^e(-/3*,^*). 



where 



(2.12) 



One can observe that A: oscillates when /? G (— /3*,/3*); contrary to the case /3 < — /3*. The function 
K may give us an intuition of the properties of the solutions of (1.7), and is useful in determining the 
behavior of the function d (see (3.6)) near the boundary of its domain. 

Theorem 2.4 There exist no solutions in H'^(M.) of equation (1.7) if any of the following conditions 
hold. 

ft) c'>l and ^k '^^'^-^'I^'^^''- -^. 

(a) F{u) > for all u, c^ > I and l3>0. 

Proof. Suppose ip e if^(M) is a solution of (1.7). Multiplying the equation by xip' and integrating 
yields the Pohozaev identity 



3((^")' - I3{^'f - (1 - c2)(^2 ^ 2F{ip) dx =. 0. 
The identity I{ip) = K{ip) may be written 



(2.13) 



(2.14) 
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Together these give 

(3p + 5) / {ip"f dx - (p + 3)/3 / {ip'f dx - (p - 1)(1 - c^) I ip^ dx^ 0. 



The term on the left side of this equation will be positive, a contradiction, when condition (i) is 
satisfied. Next, eliminating the ip" terms in the equations above gives 



2p / (v?') da; - 4(1 - c^) / ip^ dx = -(3p + 5) / F{ip) dx. 

JR JR Jr 

The conditions in (ii) imply that the left hand side is non-negative and the right hand side is negative. 

n 



3 Stability 

In this section we establish that the set of ground state solitary waves is stable under a suitable 
convexity condition. 

Theorem 3.1 (Local Existence) Suppose p > 2. Let uq — {uq,vq) G X , then there exists T > 
and the unique solution it — {u,v) e C([0,T); ^) of (1.6) such that u{0) — uq- Moreover u satisfies 
E{u) = E{uo), Q{u) = Q(uo), Qi{u 

E{u) = E{u,v 

Q{u) = Q{u,v 
Qi{u) = Qi{u,v 
Q2{u) = Q2{u,v 
Qsiu) = Q3{u,v 



h(^*o), Q2iu) = Q2{uo) and Qsiu) = Qsiuo) where 

(u^ - l3ul + u\^ + v^) - F{u) dx, 



2 

uv dx, 

u dx, 
V dx, 
udl^v dx. 



ken. 



(3.1) 

(3.2) 
(3.3) 
(3.4) 
(3.5) 



and F' ^ f and F{Q) ~ 0. Furthermore T = -|-cxd, or T < +oo and 



lim IImII^ ~ +00. 

t->T- 



Proof. First write the system (1.6) as 



wt = Bw + g{w). 



where 



B = 







dx 



dx + pdl + dl 



g{w) = {{),- f{u)x). 



The result then follows by classical semi-group theory [20, 22], once we show that B is the infinitesimal 
generator of a Co-semigroup of unitary operators on 2^ , and that g is locally Lipschitz on ^' . Define 
an inner product (•, ■) „ on ^ by 



{{ui,Vi),{u2,V2)) p = i {Ul)xx{u2)xx- l3(ui)x{u2)x+UiU2+ViV2dx. 



Then for and w ~ (u, v) £ ^ , we have 



{BvS, W)a = / V,^,^xUx,j, - iSVxxUx + VxU + {Ux + (iUxxx + Uxxxxx)v Ax 
JR 

= 

and therefore B is skew adjoint with respect to this inner product. It then fohows from Stone's 
Theorem that B is the infinitesimal generator of a Co-semigroup of unitary operators on JT. Now let 
wi,W2 G ■^ ■ Then 

\\g{w2) - 9{wi)\\ .r = \\[f{ui)- f{u2)]x\\mm) 

= \\f{ui){ui)x - .f'{u2){u2)x\\L^R) 

< ll/'("l)('"l " U2)x\\l^R) + \\{u2)x[.f'{ui) - /'(w2)]||l2(R) 

To bound the first term, we use the homogeneity of / and the imbedding of i/^(M) into L°°(IR) to 
obtain 

ii/'K)iil~(«) < ciiuiii^-\«) < cii«iii^-.;„) < c\\w,r^\ 

and thus 

||/'(ui)(wi - U2)x\\l-^ < C\\wi\f^^\\{ui - U2)x\\l2{R) < C\\wi\\^^^ \\wi - W2\\X- 

For the second term, we again use the homogeneity of / and the imbedding H^{M.) into i°°(M) to find 

\\{u2)x[f'{ui) -/'(u2)]||i2(R) < C||u2||//2(js.)(||ui||^2(R) + ||u2||h2(r))P"2||wi --U2||l2(r) 

< C'|lw2||5r(||wi||.r + \\w2\\xY~'^\\wi -W2\\.%- 

Hence g is locally Lipschitz on ^, and the proof of local existence is complete. The conservation laws 
then follow by differentiating each quantity with respect to t and using the system (1.6). D 

Definition 3.1 We say that a subset S ^ ^ is .^-stable if for every e > there exists some S > 
such that whenever 

inf |||wo ~if\\,s: -.if e S\ < S, 

the solution w of the system (1.6) with w{0) — wq exists for all t > and satisfies 



supinf |||w(i)-V!'|l^ : -0 e S*} < e. 
t>o ^ J 



Otherwise we say the set S is ^ -unstable. 

In this section we show that the stability of the set of ground states is determined by the convexity 
of the function 

d{c) = E{0) + cQ{0) (3.6) 

where — {(f, —cip) and ip G $^(/3, c). 

Tiieorem 3.2 Denote ^(/3, c) = {tp — {ip,~ap) : ip £ ^(/3,c)}. Suppose (? < \ and /3 < /3* = 
2Vl-c2. Ifd"{c) > then ^{(3, c) is ^-stable. 

Before proving Theorem 3.2, we state the basic properties of the function d. We first note that, for 
any w = {u,v) £ S^ we have 

E{w) + cQ(w) = -/(u) 5_i^(y) + 1 / (cm + vf dx. (3.7) 

2 P + 1 2 Jr 



Applying this to ip — {ip, —ap) where ip g ^(/3, c) and using the fact that /((p) = K{ip), we have 

E{0) + cQ{p>) = \l{ip) - :^K{^) = ^§^^M- 

By relation (2.4) this implies that 

d(c) = ^|^™(/3,c)?^ (3.8) 

so d is well-defined, and the properties of d may be deduced by studying the properties of the function 
7Ti(/3, c). By reasoning similar to that in [15] we obtain the following. 



Lemma 3.1 On the domain D = {(/3, c) : c^ < l,/3 < 2\/l — c^}, d is continuous and strictly 
decreasing in hath \c\ and (3. For each fixed /3, dc{P,c) exists for all but countably many c, and for 
fixed c, dpiyP^c) exists for all hut countably many j3. At points of differentiability we have 



dc{/3,c) = Q{ip) = -c (p dx 

dp{P,c) = -2 I vl da; 

for any p £ 'i^{f3,c). 

For the remainder of this section we fix /3 < 2 and regard d as a function of c only. We denote by 

Ue = Ub,c-(. = < w e ^ I inf ||?Z; - p\\x < e 



the e-neighborhood of the set of ground states 

Lemma 3.2 For each {f3,c) G £>+ = {(/3,c) : < c < l,/3 < 2\/l — c^}, there exists e > such that 
the mapping c : C/^ — > M defined by 



c{w)^c(u,v)^d ^ ( 2(„^,i) -^(")) 



is continuous. 



Proof. Since d is monotone decreasing and continuous, it follows that for fixed /3 < 2 its inverse with 
respect to d, d~^, is defined and continuous in some i5-neighborhood of d{c). It therefore remains to 
show that ^^~ri\K{u) lies in this neighborhood when u d U^ and e is sufficiently small. First observe 

that for any ui,U2 E H'^{M.) we have 



\Kiui) - Kiu2)\ <ip+l) \F{ui) - F{u2)\ dx 

JR 

= {p+l) / \f{fi{x)ui + {1- fj,{x))u2)\\ui-U2\dx 
= (p+l) / C\n{x)ui + (1 - fi{x))u2)\P\ui ~U2\dx 



< C'(II'"i|Up+i + II"2||lp+i)''II"1 - "2||lp+1- 

Thus by the embedding of H'^{R) into LP+\R), it follows that K is locally Lipschitz on H'^{R). Gr 
any p S ^(/3, c) the coercivity condition (2.3) and relations (3.8) and (2.4) imply that 

II^IIh^(m) < c-'i{p) = c-i^^^^d(c). 



Hence the set of ground states ^(/3,c) is a bounded subset of S^ . Consequently the neighborhood 
C/e is bounded for any e > 0. Thus since 2^JqriT^(v) — d{c) for any ip S ^{(3,c), the Lipschitz 
continuity of K and boundedness of U^ imply that 2('~+i) K{u) lies in the ^-neighborhood of d{c) for 
all w = {u,v) E Ue a e > is small enough. D 

Lemma 3.3 Suppose d"{c) > 0. Then there exists some Ec > such that for any 
ff e ^(/3, c) and any w G U^ we have 

E{w) - E{0) + c{w){Q{w) ~ Q{0)) > ^d"{c){c{w) - cf . 

Proof. Using Taylor's Theorem and the fact that d'{c) — Q{0) we have 

d(ci) = d{c) + Q((^)(ci - c) + irf"(c)(ci - c)2 + o(|ci - c\f 

for ci near c. Thus for ci is some (5-neighborhood of c we have 

rf(ci) > d(c) + Q{0){ci -c) + ^d"(c)(ci - c)2. 

By Lemma 3.2 it then follows that for sufficiently small tc and w = (u, v) G [/^^ we have 

d{c{w)) > die) + Qi0)ic{w) -c) + ^d"{c)ic{w) - cf 

= E{0) + cQ{0) + Qm(c{w) -c) + \d"{c){c{w) - cf 

= E{0) + c{w)Q{0) + \d"{c){c{w) - cf. 

Next suppose ip e ^(/3,c(w)). Then I{ip; l3,c{w)) — K{ip) — _^ ' d{c{w)) = K{u) and -0 minimizes 
I{-;f3,c{w)) subject to the constraint K{-) = K{u). By (3.7) we have 

1 1 If 

E{w) + c{w)Q{w) = -I{u; P, c{w)) K{u) + - {cu + vf dx 

2 p+l 2 Jr 

>^I{^;l3,c{w))-^K{ij) 
= didw)). 
Combining these inequalities proves the desired result. D 

Proof of Theorem 3.2. Suppose '^{P,c) is JT-unstable, and choose initial data Wq such that 

This implies that there exist (/3fc G ^(/3, c) such that 

lim ||u;S-¥?fclk=0. (3.9) 



k- 



'OO 



Denote by w^{t) the solutions of (1.6) with iZ;'^'(0) — Wq. Then there exist some S > Q and times tk 
(for each fc > j) such that 

inf \\w''{tk)-ip\\s: ^d. 

10 



Without loss of generality we may also suppose that S < ec and therefore w^^tk) e U^^, so that Lemma 
3.3 implies 

Eiw^t"^)) - E{0) + c{w\tumQ{w\tk)) - Qm > ^d"icMw\tk)) - c]^. (3.10) 

Next, using the fact that E and Q are continuous on <^ and conserved for solutions of equation (1.6), 
we have from equation (3.9) that 

lim Eiw'^itk)) - E{0k) = lim E{w^) - E[0k) = (3.11) 

fc— foo fc— foo 

and 

lim Q(w;''Xifc)) - Q{^k) = lim Q(4') - Q{^k) = 0. (3.12) 

fc— foo /c— >-oo 

By Lemma 3.2, the sequence of scalars c{w''{tk))) is bounded, and thus equation (3.10) implies that 

lim c{w (tk)) = c. 

By continuity of d, this implies that K{u^{t}S)) = _]^ ' d{c{w^ {tk))) converges to Jj^ d{c) . Using 
the relation (3.7) and the fact that d{c) ~ E{0k) + cQ{0k), it follows that 

ll{u''{tk)) = E{w\tk)) + cQ{w\tk)) + ^X(w^ifc)) - \ [Ricu^tk) + vHtk)f 
2 p+l 2 J 

< E{w\tk)) - Eiifik) + ciQiw'^itk)) - Q{0k)) + -^K{uHh)) + d{c) 

p+l 

-> -d{c). 

P- 1 



Thus 



.ku ^^^2(p+l) 



limsup/(u''(tfc)) < -^ -d(c), 

fc^oo P - 1 



which implies that u^{tk) is a minimizing sequence for A = _]^ d{c). By Theorem 2.1, there is a 
translated subsequence, renamed u''{tk), that converges in iJ^(IR) to some ip G ^(/3, c). To control the 
second component of w'^{tk) observe that 

1 [Ricu\h) + v\tk)f = -lliu\tk)) + -^K{u\tk)) + E{w\tk)) + cQ{w\tk)) 

2 J 2 p+l 

-> -d{c) + d{c) = 0. 
Hence v^{t].) converges in Li^iM) to —c^p^ and thus w'^{tk) converges in ^ to (p = {(f, —cip). Therefore 

inf \\w''itk)-0\\.T = O, 
a contradiction. This completes the proof of Theorem 3.2. D 



4 Instability 

In this section we establish conditions that imply orbital instability of solitary waves. 
The following theorem is a key point in the proof of the instability. 

Theorem 4.1 Let A'^uq and A^wq be in L^{M.). Assume that \f{s)\ = 0{\s\p) and \f'{s)\ = 0{\s\P''^), 
as s — >■ oo, for p > I. Suppose also that u~ {u,v) is a solutions of (1.6) with u{0) = uq. Then 
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(^) ^fp>2, 



(li) ifl<p<2, 



sup 



u{z, t) dz 



<c(l+t2/3 + t4/5) 



sup 



u{z^ t) dz 



<^Al-(p-l)/3^^1-(p-l)/5\ 



for t G (0, T), where T is the maximum existence time for u, and the constant C > depends only on 
\\y.o\\x, I awdsuptg[o,T) ^(i)!!^^- 

To prove Theorem 4.1, a series of useful lemmas are laid out. The first one is the well-known Van 
der Corput lemma [23] as follows. 

Lemma 4.1 Let h he either convex or concave on [a, b] with — oo < a <h < +oo. Then 



Jhii) 



d? 



<A{mm\h"{0\ 



-1/2 



ifh" ^0 in [a,b]. 

Lemma 4.2 Suppose h is twice differentiable on R and 

(i) h" has finitely many zeroes, all of which are of order qi or less. 

(ii) there exist positive constants Ci and C2 such that \h"{^)\ > C2|^|^^ whenever \^\ > Ci. 
Then there exists a constant C such that 



for < t < 1, and 



fort > 1. 



,ith(c) 



dC 



< Ct-^/i^+11) 



Proof. First suppose < i < 1. Given R > Ci, set f7i = {^ : |^| < R} and f^s = {? : |^| > R}- 
Then \h"{^)\ > C2R'''' for ^ e f^a, so by Lemma 4.1 we have 



ei*M?) d^ 



< CR-'i^'H-'^/^, 



while on J7i we have 



,ith(c) 



de 



<R. 



For < t < 1 sufficiently small, we may set R = t^ 1/(2+92) qj^^ ^j^g result follows. 

Next suppose t > land let ^1, . . . , ^n denote the zeroes of h" . For e > let /fe = {^ : |^ — ^a:| < e} 
for each fc, and set fii = |J^ Ik and VI2 = M\r2ii. Then we have 



pi*''(«) 



dC 



<Ce. 



Since each zero of h" is at most order gi, there exists C3 > such that for e > sufficiently small, we 
have \h"(C)\ > Cje'^^ for £_ e ^2- It then follows from Lemma 4.1 that 



Jthii) 



dC 



< ct-^/h-'''/^ 



For t sufficiently large, we may set e = i 1/(2+91) ^^d the estimate follows. 
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Lemma 4.3 Let P <2 and set h{^, a) = ^^^ „ ^^4 ^ ^6 ^ ^^^ 
(i) If /3 ^ 0, there exists a positive constant C such that 



sup 



gitft(«,") d^ 



< Ci-1/3 



for all t > 0. 

(ii) If p — Q there exists a positive constant C such that 



sup 

agR 



ei*''(C.") d^ 



<c{t-''^ + t 



1/3 , ^-l/5^ 



for all i > 0. 
Proof. First observe that h is an even C°°-function in M \ {0} with 

d'^h _ Id (-3/? + (2/32 + 10)^2 - 9/?e* + 6C^) 



(4.1) 



a^2 (1_ ^^2 +^4)3/2 

Since the polynomial -3/3 + (2/3^ + 10)^^ - 9/3,^* + 6^'' is increasing in (^ for /3 < 2, it follows that 

(i) if /3 > 0, /i" has three simple zeroes, 0, ^0 > and —^0, 

(ii) if /3 < 0, h" has one simple zero, ■^ = 0, 

(iii) if /3 = 0, h" has a zero of order 3 at ^ = 0. 

In cases (i) and (ii) the result then follows from Lemma 4.2 with qi — q2 = 1, while for /3 = it follows 
from the same lemma with gi = 3 and 52 = 1- D 

Lemma 4.4 If u e LP{M), 1 < p < 00, then A^^u G Lp{R) and ||A^2u||ip(R) < C||u||iP(K), for some 
C>0. 

Proof. The proof follows from Young's inequality and the fact G{x) = exp(— |a;|) £ L^(M), where 
A-^u = G*u and G{0 = {I + ^^)~'^. □ 

The following lemma gives a time estimate on the solutions of the linearized problem. 
Lemma 4.5 Let S{t) be the Cq group of unitary operators for the linearized problem of (1.6) 



-1 

-1 - f5dl - dl 



Ut+ \ 1 «a2 n4 A ) ^2; — 0, 



with u(0) = Uo(Mo,fo)- //A^wo G £^(M) and v^ e ^^(M), t/ien S{t)uQ G L°°(M) x L°°(M) and 
||5'(i)uo||L~(R)xL~(R) < C (i-1/3 + i-1/^) (IIA^uoIIlmr) + II^'oIIli(R)) : 

where C > is a constant. 

Proof. Since 

cosit^iSiO) 5^sin(tC^(0) 

71(f) = Sit)uo{x) = / e'-« I I MO de, 

ii?(Osin(ie^(0) cos(tC^(0) 
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where '§{$,) = -^/l — /^^^ _|_ ^4^ j^ jg deduced from Fubini's theorem and Lemmas 4.3 and 4.4 that 



uo ± A-2„o ) e"^^''^^)*-/*) de 



2{t)\ = \S{t)M^)\<Y^ f 

JR 

<V f \uo{z)±A-'vo{z)\ f e' 

Js. Jm 

+ E [ \voiz) ± A-\oiz)\ [ e' 



^±^tro)e"«w«)±-/*)de 



E 



(.ro±^(0^)e"«(*(«)±^/*)d^ 



E 

t5(lJ(0±:r/t-2/t) j^ 
tC(.9(0±x/t-z/t) 



t;o ± A-2wo ) e'**^''^^)*-/*) d^ 



dC 



dx, 



where the sums are over all two sign combinations. Therefore, we obtain from Lemma 4.3 that 



\u{t)\ < sup 



Jth{i.a) 



d^ 



(ll'"ollLi(R)xiHK) + 11^ ^'oIIli(K) + 11^ ■"oIIli(R)) 



for t > 0. Hence, for some C > 0, it is concluded 

\uit)\ < c (t-i/^ + t-i/s) (II^oIUhk) + II A'^oIUh 
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A proof of Theorem 4.1 is now in sight. 
Proof of Theorem 4.1. Let w{t) = S(t)uQ, then w satisfies 



Wt 



-l-(idi-dt 
Then the solution u{t) of (1.6) can be written 



2 q4 n ) ^^ = 0' "^ith w{Q) = Mo- 



We should estimate 



where 



u{t) = w;(i) - 9, y^ ^(t - r) (^ ^^^^^^^ 1 dr. 



^w = ^w-X^(^--)Uk-)) ""' 



U{x,t) = 
First observe using (4.2) that 



u{z,t) dz and iy(a;,t) 







w{z^ t) Az. 



(4.2) 



(4.3) 



where C/o(x) — J_ uo{z) dz. Now Lemma 4.5 implies that 

\W{x,t)\<\\uo\\mM)xLHR) + {\\A^ua\\mji) + \\A\o\\mR)) / (r-i/^ + r-i/^) dr 



< 



(IIA^uoIUme) + \\A\o\\mm) (i + 1'/' + 1^^' 
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Setting 



' 



F(x,t)= / S{t-T)\ ^^^^^^^ ) dr, 



and using Lemma 4.5 again, it follows that 



y{^M < f ((i - r)-"' + (t - r)-'") \\f{n{T))\\L^m dr. (4.4) 



/o 
On the other hand, it is deduced from Cauchy-Schwarz inequality that 

\Y{x,t)\< f f (l+ , l^^ =1 |/(^^))|dedT< r||/W-,^))lk=(M)dr. (4.5) 

Jo Jr \ V^~ P^ +i J ■'0 

Since H^{R) ^ L°°{R) and |/(s)| = 0(|s|p) and |/'(.s) = 0{\s\p-^)\ as s ^ 0, for p > 1, it transpires 
that ||/(m)||li(r) ^ C*; provided p > 2, for some positive constant C which depends only on / and 
suPfg[o,T) \\^{t)\\.r- Hence, if p > 2, 

\Yi^, i)\< f [{t - r)-'^' + it- r)-'/') dr < C (t^/^ + ^4/5^ _ 

If 1 < p < 2, it is straightforward to check that ||/(u)||^i,2/p(r) < C, for some C > 0. Since (4.4) 
and (4.5) hold for any / S L^(E) n i/^(E), a straightforward interpolation thus can be applied for the 
mapping S{t — r) as in (4.4) and (4.5). Thus the same argument proves that 



,p-i 



\Y{x,t)\<J^ ((t-r)-V3 + (i_^)-i/5y- ||/(u(,r))||^..v.(M)dr 

<^/^l-(p-l)/3^^1-(p-l)/5V 

By combining the estimates of Y and W, the proof of Theorem 4.1 is now completed. D 

Given (f £ '^{f3,c) and e > 0, we define the "tube" 

f^<3.e = i w € JT; inf \\v -u\\x <e. 

and the operator 

H = L"{0) = E"i0) + cQ"i0). 

The main instability result is the following. 

Theorem 4.2 Suppose c^ < 1 and (3 < (3^,. If there exists tp £ ^ such that d^Tp G L-^(M) x L^(M) and 
^e H'^iR) X H^{R) 

1. (Q'(^),dj)^G, 

2. {HdJ, dj) < 0, 

then O0 is ^Sf -unstable. 

Lemma 4.6 Let c^ < 1 and P < j3i, and ip G ^(/?, c) be fixed. There exist eo > and a unique C^ 
map a : ^ip,to — > K such that a{(p) — 0, and for all it G ^ip,eo o,nd any r G M, 

(i) {u{- - a{u)), d^ip) = 0, 
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(ii) a{u{- + r)) = a{u) — r, 
(ill) a'{u) = (^u,dl0(^-a(u))) ^^'?i- - "(^))' ""'^ 
(iv) {a'{u),u) = 0, if u ^ O^. 

Proof. The proof follows the line of reasoning laid down in Theorem 3.1 in [12] and Lemma 3.8 in 

[19]. D 

Let ip be as in Theorem 4.2. Define another vector field Br by 



Bj;{u) = Jfdxijj{- - a{u)) 



(u,d.j,'4){- - a{u)) 
{u,dlLp{--a{u))) 



Jtdlip{--a{u)), 



for u e ^if,ii 7 where J(f = _. n ] ■ Geometrically, Br can be interpreted as the derivative of the 
orthogonal component of t_^/.\ijj with regard to T_^i.\d,j.0. 

Lemma 4.7 Let ip be as in Theorem J^.2. Then the map Br : r^i^e^ —^ 3^ is C^ with bounded 
derivative. Moreover, 

(i) Br commutes with translations, 
(i^) {B^{u),J(fu) =0, ifue ^0^eo, 
(Hi) B^{0) = dx.J^ip, if {0,dxi^} = 0. 
Proof. The proof follows the same lines from the proof of Lemma 3.5 in [1] or Lemma 3.3 in [2]. D 



Before starting with the proof of Theorem 4.2, we state and prove the following lemma which shows 
the boundedness of the Liapunov function (see (4.13)). 

Lemma 4.8 Let tp be as in Theorem J^.2, ug = ("Oi'^^o) ^e in Hi^gj such that Auq, Awq G L'^(]R) and f 
satisfy the assumptions of Theorem 4-1- If u{t) is a solution of (1.6) which corresponds to the initial 
data uo and u{t) G ^f,e3, for t e [0,T], then 



'il'{x — a{il{t))) ■ u{x, t) dx 



<c(l + t2/3+t4/5) 



(4.6) 



for t e [0,r), where T is the maximum existence time for u, and the constant C > depends on 
||A^uo||li(R), ||A^uo||li(r), / andLp. 

Proof. Let H be the Heaviside function and 7 = /jg dxtp{x) dx. Then the left hand side of (4.6) may 
be written 

I t/i (x — a{u{t))) — 7ElI(a; — a(u{t))) ) • u{x, t) dx + •y ■ / u{x, t) dx. 

So it follows from Cauchy-Schwarz inequality and Theorem 4.1 that 

'\\}(x — a{il{t))) ■ u{x, t) dx 



< 



ip-jM 



L2(R)xL2( 



\um 



L2(M)xL2( 



+ C (l + ^2/3 + ^4/5^ 
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We show that tp — 7H e L^(M) x L^(M). Indeed, Minkowski's inequahty yields that 



ip-'jM. 



< 



L2(K)xL2(R) 



.1/2 

\rP{x)\'^dx] + 



+00 \ 1/2 

\^{x) ^ ^m{x)\^ dx 




2 



+00 



dx-)p{z) dz 



1/2 



dx' 



< ||V'I|l2(k)xL2(r) + / \dxip{x)\\/\x\ dx 

Jr 

< \\'>P\\h^(R}xH^{R) < +00. 



Hence, for t e [0,T), we obtain 

iplx — a(u{t))) ■ u{x, t) dx 



<c(l + t2/3+t4/5y 
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All the elements are now in place to prove the instability result in Theorem 4.2. 
Proof of Theorem 4.2. First we claim that there exist 63 > and (73 > such that for each 

L{0) < L{uo) + ^{uo)s, (4.7) 

for some s G (— CT3,o-3), where ^{u) — {L'{u),Br{u)). 

For uq £ ^0,eo, where eg is given in Lemma 4.6, consider the initial value problem 



£u{s) = B--{u{s)) 
u(0) = uo- 



(4.8) 



By Lemma 4.7, we have that (4.8) admits for each uq G ^i^.^o ^ unique maximal solution u G 
C^((— cr, o');r2^^e„), where a G (0, +00]. Moreover for each ei < eo, there exists cti > such that 
""(wq) ^ 1^1: for all Uq G ^i^,ei- Hcncc we can define for fixed ei, cti, the following dynamical system 



^ V,co 



-^ : (-0-1, cr) X f^^^ei 

(s,uo) i~> '^(s)wo, 

where s — >■ '^{s)uo is the maximal solution of (4.8) with initial data uq. It is also clear from Lemma 

4.7 that '^ is a C"'^— function, also we have that for each uq G ^ip,ei, the function s — > ^{s)uo is C^ 

for s G (— cri,cri), and the flow s — > '^(s)'Uo commutes with translations. One can also observe from 

the relation 

pt pt 

P{s)Ta(''if(s)0)dlifds 



'^{t)tp = (p+ Ta(^ir(s)'^)dx1pds 



that ^{s)(pe i7'^(M), r > 3/2, for ah s G {-ai,ai), where 

^(s)(^,T„(^(t)^)a^v^ 



p(^) 



(^(t)^. 



^C^Wv 



?)5^<^) 



Now we obtain from Taylor's theorem that there is p G (0, 1) such that 

L(^(s)uo) = L{uo) + ^{uq)s + -R{^{qs)uq)s', 
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where R{u) = {L"{u)B^,B^{u)) + {L" {u),B'Ju){B^{u))). Since R and ^ are continuous, L'{Lp) == 

and R{0) < 0, there exists £2 G (0, ei] and (72 G (0, Ui] such that (4.7) holds for uq £ B{0,e2) and 
s £ (— fT2,o'2)- On the other hand, it is straightforward to verify that 



P{'^is)uo] 



= and -— P('^(s)wo] 

{uQ.s) = {v,0) d* 



= (p'(v?),a,V), 

(iTo,s) = ('P,0) 



where P is defined in Theorem 2.2. We show that {P'{(f), dxi^) 7^ 0. Otherwise, 93,-0 would be tangent 
to ,y\^ at 0, where ,yV is defined in Theorem 2.2. Hence, {L"{(f)dxtjj,dxtp) > 0, since (p minimizes 
L on ^ by Theorem 2.2. But this contradicts Theorem 2.2. Therefore, by the implicit function 
theorem, there exist £3 € (0,62) and (T3 € (0,172) such that for all uq € B(p, 63, there exists a unique 
s = s{uo) G (—0-3, a^) such that P{'^{s)uo) = 0. Then applying (4.7) to {uq, s{uo)) £ B(p, e3X (— 0-3, 0-3) 
and using the fact (p minimizes L on ^ , we have that for uq £ B0, 63 there exists s £ (—173, 173) such 
that S{0) < L{^{s)uo) < L{uo) + ■^{uq)s. This inequality can be extended to ^(p^f,^ from the gauge 
invariance. 

Since '^{s)uq commutes with r,,, it follows by replacing uq with '2^(s)'Uo in (4.7) and then 5 — —s 
that 

L{0) < L{'^{5)0) - ^{'^{6)0)5, (4.9) 

for all 5 £ (—(73,(73). Moreover, using Taylor's theorem again and the fact [^{0) = 0, it follows that 
the map & 1— > L{'^{5)0) has a strict local maximum at (5 = 0. Hence, we obtain 

L{^{5)0)<L{0), 5^0, (5 €(-(74, (74), (4.10) 

where CT4 £ (0,0-3]. Thus it follows from (4.9) that 

^('^^((5)^) <0, 5 £(0,(74). (4.11) 

Now let 5j £ (0,(74) be such that 5j — > as j — ?> 00, and consider the sequences of initial data 
uoj ~ ^{^j)0- It is clear that Uqj £ H^{M.), r > 3/2 for all positive integers j and uqj -^ (p in J^ 
as j — > 00. We need only verify that the solution Uj{t) — ^{t)uQj of (1.6) with Uj{0) — uqj escapes 
from 51^ g., , for all positive integers j in finite time. Define 

Tj = sup{t' > 0; Uj(t) £ f7^,,3, Vi £ {0,t')} 

and 

^ = {u£n0^,,]L{u)<L{ip), ^{u) <0}. 

It follows from (4.7) that for all j £ N and t £ {0,Tj), there exists s — Sj(t) £ (—(73,03) satisfying 
L{ip) < L(moj) + ^{uj{t))s. By (4.10) and (4.11), mqj £ 9\ and therefore Uj(t) £ ^ for all t £ [0, T^]. 
Indeed, if ^{uj{to)) > for some to £ [0,Tj], then the continuity of 3^ imphes that there exists some 
ii £ [0,Tj] satisfying 3^{uj{ti)) — 0, and consequently L{(p) < L{uoj), which contradicts uo.j € ^• 
Hence, & is bounded away from zero and 

_^(^r^.)>iiM^iii^^,^. >o, vie[o,r,]. (4.12) 

Now suppose that for some j, Tj — +(X). Then we define a Liapunov function 

A{t) = / 0(x - a{uj)) ■ Uj{x,t) dx, t £ [0,Tj]. (4.13) 

Since 

^ = dxJTE'iu,), 
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^ = /a'{u,{t)), ^\ (dj{- - a{u,{t))), u,{t)) + U{u,{t))), ^ 



then we have 



= - {B^{uj{t)),L'iuj{t)))+c{B^,{u,it)),Q'{uj{t))) - -^(u,(i)), 
for t e [0,Tj]. Therefore it is deduced from (4.12) that 

-^>r?, >0, Vie[0,T,]. 

This contradicts the boundedness of A{t) in Lemma 4.8. Consequently Tj < +00 for all j, which 
means that Uj eventually leaves ^^1^,63. This completes the proof. D 

The remaining results of this section are applications of Theorem 4.2. In verifying the hypotheses 
of this theorem, we will use the fact that for any wi = (ui, vi) and W2 — {u2, V2) in ^ we have 

{Hwi,W2) = / «" + M' + (1 - c^)"i - f{^)ui)u2 + (cui + vi){cu2 + V2) dx . (4.14) 

In view of this, we define Hi = d^ + f3d'^ + (1 — c^) + f'if)- Our first result is the following complement 
of Theorem 3.2. 

Theorem 4.3 Suppose (3 < 2 and assume there exists a C^ map c ^-^ ip ^ ^(/3, c) for c < Ci,. If 
d" (c) < 0, then O^ is ^ -unstable for any ip € ^(/3, c). 

Proof. Define 

i^ix) = / ip{y) - -JiT^'fciy) dy, 

where (pc = j;:0 — {p>c, —cfc — p)- We need to show that ij) satisfies the hypotheses of Theorem 4.2. 
Now 

{Q'{0),dj) = {Q'ip), p) - ^^Q(¥?) = 2d'(c) = 2d'{c) = 
so the first hypothesis is satisfied. To show that the second hypothesis is satisfied, first note that 

[HdJ, dj) = {Hp, p) - ^ {Hp, (p),) + (^^) {H0C, Pc) ■ 

Using the homogeneity of / and the solitary wave equation, we have 

Hiip)^fip)^fip)p^il-p)fip), 
so by relation (4.14) it follows that 

{Hip, p) = {l-p) [ f{p)p dx - (1 - p){p + 1) / Fip) dx = -2(p + l)d(c) 

Jr Jr 

and 



{Hp, p,) - (1 - p) / f{p)pc dx^{l-p)[ F{p) dx = -2rf'(c). 

JR \JR J c 

By differentiating the solitary wave equation with respect to c, it follows that 

Hipc = 2cp3, 
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so 



.^ , d'ic) 



c 
d'{c)\ d'{c) 



c 
= -d"(c). 

It now follows that 

{HdJ,dj) = -2(p+ l)d(c) - ^(-2d'(c)) + (^^y (-d"(c)) 
= -2(p+lMc) + l^<0 
since d"{c) < 0. This completes the proof. D 

We next apply Theorem 4.2 to obtain conditions on p, (3 and c that imply orbital instability. For 
our choices of unstable direction we will use the following. 

(i) ipx = (<y5, +C'/') ^ for small c, and any p > 1. 
(ii) ^x = '0 + '2xipx ~ for large p. 
Lemma 4.9 Lei iSa;-^ — + 2xipx — if + '^x'-Px, ~cif + ^xip')). Then (Q' [0), dxip) — and 

(Hdx^,dxA - (l-P)(p-3) j^(^) + /• 24(^")2 _ 4^(^.)2 d^^ 
\ / P + 1 Jm 

Proof. First, we have 

(q'{0), dxTp) = / -2cip{ip + 2cip') da; = 

as claimed. Next we have 

(dxH^^,dxA= /[(¥^ + 2V)"" + /3(¥' + 2V)" + (l-c')((/3 + 2V)-/'(^)(^ + 2V)](^ + 2V)da;, 

which may be split into three terms: 

Ai^ J w" + pf" + (1 - c2)^ - n^)^]^ dx, 

A2^2 f y" + /3(^" + (1 - c")^ - f\ip)^]{2x^') dx, 

A3= f [{2x^T" + /3{2xpT + (1 - c^){2x^') - f{p){2xp')]{2xp') dx. 

Jr 

Since ip"" + ^ip" + (1 - c^)ip - f{ip)ip = (1 - p)/((/3) we have 

A, = {\-p){p+l) [ F{ip) dx - {l-p)K{ip) 
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Figure 2: The regions of instability guaranteed by Theorem 4.4. The regions described in part (i) he 
between the upper and lower curves on the first plot. The regions described in part (ii) lie to the left 
of the curves in the second plot. Both regions grow to fill the domain of d as p increases. 



and 



A2 = -4(1 - p) I F{^) dx = -^^-^K{^). 

p+ 1 



For A^ first observe that by differentiating (1.7) we obtain if'"" + jSip"' + (1 — c^)¥'' — .f'{'p)'fi', and thus 
{2xip'y"' + (i{2xf')" + (1 - c^)(2xf') - f{f){2xf') = 8f"" + 4/3^". 

Thus 

A3= f (8^"" + 'i/3(p")2xip' dx= [ 24(^")^ - 4/3(93')^ dx 
Jr Jr 

so summing Ai , A2 and A^ yields the result of the lemma. D 



Theorem 4.4 Suppose c^ < 1, /3 < /3* = 2\/l - c^ and tp e ^(/3,c). Recall that c* = W 1 - /3+/4. 
Then 0{ip) is ^ -unstable in the following cases. 



(i) p>l andc^ < (2_i) cl. 
(n) p>9,c^<landl3< (^E§^) /3. 



Proof. To prove the first statement, consider the choice dxi/j — {(p,c(p). It is easy to see that 
Q'{(f),dxi^) — 0. Next we compute 



d^H^P, d^^p) = / [if"" + 13 f" + ip- f'{ip)f)ip + (cipf + 2apap dx 
Jr 

= [ {ip"" + /?(^" + (1 - c2)<^ - f'iip)ip)f + 4cV' dx 
Jr 

= {l-p)K{ip) + Ac^ / ip"^ dx. 



First suppose /3 < 0, in which case c* = 1. Then I{ip) > (1 — c^) /^ p"^ dx, so 



5,iJ^,9,V)< {1-P + 



4c^ 
l-c2 



^M- 
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Now suppose /3 > 0. Then 



m = / i^'r ' PWY + ^^' dx+(i-c'- ^) I ^^ dx 



> ( 1 _ C2 _ ^ ) / ^2 ^X 



4 '^ V 4 

and thus 

4c- 



Hence for any /3 < 2 we have 



and this quantity is negative when condition (i) is satisfied. 

To prove (ii), we use the choice of unstable direction given in Lemma 4.9. Multiplying the solitary 
wave equation by xip' and integrating yields the Pohozaev identity 

3((^")^ - /3{<fi'f - (1 - c2)(^2 _^ 2F{ip) dx = 0. (4.15) 

The identity I{y:>) — K{ip) may be written 

{ip"f - Pi^'f + (1 - c2)y,2 _ (p + i)F{ip) dx = 0. (4.16) 



/R 

Together these give 

jr\2 r,/^(',,'^2 



A{ip"Y - 2p{ip'y - (p - l)F(^) dx = 0. 
Together with the result of Lemma 4.9, this gives 

\ / P + 1 Jr 

Since 

/(^) >(/?,- /3) /(^')'dx 

it then follows that 

(...*>.*).(^.il^^fc^)A-M. 

The term in parentheses is negative when /3 satisfies condition (ii) above. D 

5 Further Properties of d. 

In this section we establish further properties of the function d. We first obtain bounds on the 



function d as c approaches ±c* = i/l — /3^/4. To obtain these bounds, we use trial functions to obtain 

bounds on the Rayleigh quotient that defines m{f5,c). To motivation the choice of trial function, we 
observe that solutions of the solitary wave equation (1.7) have tails that decay like solutions of the 
linear equation 

ifi"" + pip" + (l - c^)ip = 0. (5.1) 

22 



The fundamental solution of this equation is the function h defined by (2.10). Recalling that h is given 
explicitly by the expressions in (2.11), we see that h e i/^(M), and is thus a valid trial function provided 
K{h) > 0. The fact that K{h) > will be verified below. Since scaling has no effect on the Rayleigh 
quotient that defines m, we use the following scaled versions of h for simplicity. If — /3* < /3 < /?*, 
define 

u{x) — e^'^l^l(a;cos(a;a;) + cr sin(a;|2;|)) (5.2) 

and if /3 < — /3* define 

v(a;) = Ase-^il^l-Aie-^^l^l (5.3) 

where Ai, A2, <y and uj are defined by (2.12). 

Theorem 5.1 Suppose f{u) — \u\p~^u. Fix /3 < 2. Then 

d(;9,c) = o((c*-c)5oS)^ 

as c approaches c* . 



Proof. First consider < /? < 2. Then c* = \/l — 4/3^, and it follows that 



UJ^ yjjj2 + 0{c., -c) 

as c — >■ c*. For the trial function u given by (5.2), a direct calculation reveals that l(u) — 4(juj'^{a'^+uj'^), 
and by calculations similar to those in [16] we have 



Kiu) = / Wr' dx > -^ 

JR t;(o 



0{a) 
for small cr > 0. Thus 

liu) 



liU] / P+3 \ / P+3 \ 



if(u)2/(P+l) 
as C r Cj^ . 

Next, when /3 < we have c* = 1, and 

Ai — 0{y/c^ - c) 

A2 = y^ + 0(c*-c) 

as c — ?> c*. For the trial function v given by (5.3), another direct calculation reveals that I{v) — 
2(A2 — Ai)AiA2(A| — Af ), and by calculations similar to those in [16] we have 

''■<"' ^ oh 

for small Ai > 0, and thus 

as C r c^ . 

The result then follows by the relation between d and m. D 

Corollary 5.1 Suppose f{u) = \u\p^^u where 1 < p < 5. Fix /3 < 2. Then there exist c arbitrarily 
close to c* such that'^{j3^c) is ^'-stable. 
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I o p+3 

Proof. Since ^f"2\N > 1 when 1 < p < 5, the function (c* — c) ^(p-i) is convex and vanishes at c = c*. 
Thus by Theorem 5.1, d vanishes at c = c* and is bounded above by a convex function. Since d is 
positive, this imphes that there exist c arbitrarily close to c* such that d"{c) > 0, and the result then 
follows from Theorem 3.2. D 

Remark 5.1 The results of Theorem 5.1 and Corollary 5.1 also hold Jor the even nonlinearity f{u) — 
\u\P~^^ in the case that /3 < since the trial function v is positive for small Ai (a near 1). However, for 
< /3 < 2 the non-positivity of u only allows one to obtain the weaker estimate d{/3,c) — 0{^c^, — c) 
which does not imply convexity of d near c* . 

We next present the main scaling identity satisfied by the function d. 

Theorem 5.2 Let c^ < 1 and j3 < j3.j, = 2\J\ — (? . Then for any < r < (1 — c?)^^!"^ we have 

I 3p + 5 

d(r/3, a/I - r2(l - c^)) = r ^(p-i) d(/3, c). 



Proof. Recall that 



where 



to(/?, c) = inf ■' 



i^(u)2/(P+l) 

I{u) = I[u- I3,c)= I ul^ - Pul + (1 - c^)u^ dx 



and 



K{u) = (p + 1) / F{u) dx. 
Given any u G iJ^(M), we set v{x) — r^/^u{r~-^/^x). Then 

liv; /3, c) = / uL - r^ul + r^il- c^) dx = I{u- r/3, ^l-r^{l~c^)) 



and K{v) — r i K{u). If we then suppose v achieves the minimum m{l3,c) it follows that 

,^ , I(v:B,c) 3p+5 Ku: rB, ■v./l — r^fl — c^)) 3p+5 r- — — 

m B c) = ^ ' = r 2(p+i) ^ ' ^' ^ 1^ 11 > r ^(p+'^)m rB \/l ~ r^il — c^)^ 

^'^' ^ /f(t;)2/(p+i) ii:(u)2/(p+i) - TTU^rp,yi r yi c )). 

By supposing that u achieves the minimum m{rf3, ^J\ — r'^{\ — c^)) we obtain the reverse inequality, 
and the result then follows by the relation between d and m. D 

Remark 5.2 This scaling property implies that all values of d on any semi-ellipse /3 = k\/l — c? with 
k <2 are determined by any single value of d on that semi-ellipse. 

Setting /3 = and r^ = (1 — c^)^^ in Theorem 5.2 gives the following result. 

Corollary 5.2 When /3 = 0, d{c) = (1 - c^)^^^d{0), and it follows that 
(i) Ifp > 9, then d"{c) < for c^ < 1, 
(ii) Ifp < 9, then d"{c) < for c^ < ^^^ and d"{c) > for c^ > 2(p - l)p + 7. 

Theorem 5.3 Suppose d is twice differentiable on its domain, and consider the curve F^ = {(/3, c) : 
< c < 1, /3 = k\/l — c^} for some k < 2. Then dcdfi, c) changes sign at most once along Ffc. 
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Figure 3: Plots of d{c) for /? = and p = 2, 3, . . . , 10. The values of d{0) were found via the numerical 
methods described in Section 6. 



Proof. We present two proofs of this fact, both of which make use of the scaling property of d. First, 
setting r = (1 — c^)~^/^ in Theorem 5.2 gives 



d{/3,c)^{l~cyd{P/Vl-c^0) 
where 7 = ./_., . Equivalently, setting s — \/l — c^ we have 

d(/3, Vl-s2) = s^'^d(/3/s, 0). 
Differentiating once with respect to s gives 



or equivalently 

4(/3, v/l-s2) = Vl-s^ [-27s2T-2d(/3/s, 0) + /3s27-3^^(^/^^ 0)] 
Differentiating again with respect to s then gives 



4c (/3, Vl-S2) 



vr^^ \/r 



-27527-2^(^/5^ 0) + /3s2^-3d/3(/?/s, 0)] 



+ Vl-s2 [-27(27 - 2)s4^^3^(/?/s, 0) 

+/3(27 - Z)s^-^dp{P/s, 0) - /32s2t-54;3(/3/s, 0)] . 



25 



Now denote /3o — jij \f\^-^ — P/s. Then this becomes 
dec(/3, c) = -2-fs^^-^d{Po, 0) + /3s2T-3rf^(/3o, 0) 

+ (1 - s2) [27(27 - 2)s27-4rf(^^^ 0) ~ ^(47 - 3)s2^-5^^(^^^ 0) ^ ^2^27-6d^^(/3o, 0)] 

+ c^ [27(27 - 2)s2^-4d(/3o, 0) - /3o(47 - 3)s2^-4d^(/3o, 0) + pls^^^-^d^p^P^, 0)] 

Simphfication yields 

4c(/3, c) = s^^-^ [27d(/3o, 0)(c2(27 - 1) - 1) + MpiP^. 0)(1 - (47 - 2)c^) + c^ (ildpp{po,Q)\ 
= 5^7-4 [c2(27(27- l)d(/3o,0) - 2(27- l)/3od/3(/3o,0) + /3o'd0M/?o,O)) 
+(/3odMA,0)-27d(/?o,0))] 

Since the bracketed term is Hnear in c^, this shows that dec changes sign at most once on r^, and the 
change of sign occurs when c — vP, where 



-/3odff(/3o,0) + 27d(/3o,0) 

27(27 - l)d(/3o, 0) - 2(27 - l)l3odpWo,0) + l3^dp0il3o,O) 



^ = ^(/3o, 7) =..,.., ,,,,, ., .,.., ,,,, ,, ,, , ,2 



provided < P < 1. 

Ahernately choose any point (/?o, cq) G F^ with cq 7^ 0. Then applying Theorem 5.2 with r — /?o//? 
gives 

d{l3,c)=(-^) d{f3o,co) 

where q = 2(^~!il) ■ Differentiating twice with respect to c and using the relation cq — \/l — (3q{1 — c^)//?^ 
we have 

1 / /3 \ '"** 

rfcc(/3, c) = 3. _ 2N h^- [(1 ~ Co)coc^4c(/3o, Co) + (Cq - c^)4(/3o, Co)] . 

CqU c5j \PqJ 

The term outside the brackets is positive, while the bracketed term is linear in c^ and therefore can 
change sign at most once for < c < 1. The change of sign occurs when c = vP, where 

p = p(co) = cg4(/3o,co) 

{cl - l)codcc{(3o, Co) + 4(/3o, Co) 

provided < P < 1. D 

Remark 5.3 Theorem 5.3 does not imply that for /3 fixed dec has at m,ost one sign change as c varies. 
Indeed, when p = A there exist j3 for which dec changes sign three times as c varies from to c^,. See 
Figure 4. 



6 Numerical Results 

In this section we present numerical calculations of d and its derivatives for the nonlinearities 
f{u) = \u\^ and f{u) = \u\p~^u for several values of p. These results illustrate precisely the regions 
in the (/?, c)-plane where dec is positive and negative, hence where the solitary waves are stable or 
unstable. 
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f(u)=|u|''u, |i=1 




0.4 0.5 



Figure 4: When f{u) — \u\'^u and /3 = 1, the sign of d"(c) changes sign three times. 



The method consists of numerically computing a solitary wave ip for given (/3,c) and using the 
relations 



dW,c) 



p-1 

2{p + l 



■K{ip) 



dp{l3,c)^-- 
dcil3,c) = -c 



fx dx 



f dx 



to compute d and its first derivatives. By then doing this for several values of (/3, c) the second 
derivatives dec and dpp may be calculated numerically. By the scaling relation in Theorem 5.2, it 
suffices to perform these calculations over the segments 

5i = {(/3,c):c = 0,-l</3<2} 
52 = {(/3,c):/? = -l,0<c<l}, 



since for every k < 2 the semi-ellipse Fj. = {(/3,c) : < c < 1,(3 = ky/l — c^} passes through either 
5*1 or 82- The calculations in the proof of Theorem 5.3 may then be used to determine the locations 
where dec changes sign. 

To compute the solitary waves, the following spectral method due to Petviashvili. The Fourier 
transform of the solitary wave equation (1.7) is 



so we perform the iteration 



r - pe + (1 - c'))0 = fiv) 

^'+' ^4 _ ^^2 + (1 _ c2) 

where the stabilizing factor M is given by 



M 



S^I{Vk)'f>kd£, 



27 




Figure 5: The domain of d, {(/?, c) : c^ < 1, ^ < 2Vl — c^}. Also shown are the semi-eUipses Tk along 
which the scaling relation determines the values of d, and the segments 5*1 and S2 along which the 
numerical calculations were performed. 



The convergence properties of this method were studied in [21], where it was shown that the exponent 
^Y of the stabilizing factor M yields the fastest rate of convergence. In the case of the nonlinearity 
/(u) = uP for integer p, there exist exact solutions of the form 



ip{x) 



(p + 3)(3p + i; 
8b + 1) 



sech^-i 



p-1 

4(p+l) 



v/-(p2 + 2p + 5)//3i 



when /3 = — ( ^ — I — |-j- J \/l — c^ ([6]). On the spatial domain [—200, 200] the numerically computed 

solitary waves very closely approximate the exact solutions, with an L^ error on the order of 10~^ after 
about 100 iterations using Gaussian initial data. 

The results of these computations for the odd nonlinearity f{u) — \u\p~^u and even nonlinearity 
/(u) = |m|^ are shown in Figures 6 and 7, respectively. Each curve corresponds to a different choice of 
the power p, and separates the domain D^ into two regions 

i?„ = {(/3,c)ei5+:4c(/3,c)<0} 
D, = {(/3,c)eI?+:4c(/3,c)>0}. 

Since dcc(/3,0) < for all /?, the region of unstable solitary waves, D„, is the "lower" region that 
contains the /3-axis, while the region of stable solitary waves, Dg, is the remaining region. Several 
observations may be made regarding the stable and unstable regions. 

Observation 6.1 (i) For p < 5, the stable region Dg is unbounded and for each fixed /3 contains 
points {l3,c) near (/3,c*), in agreement with the result of Corollary 5.1. 

(ii) For p > 5, the stable region Dg is bounded, and when p > 5 appears to consist of the set of points 
interior to a smooth closed curve that passes through (0, 1). 
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Figure 6: Nodal sets of dec for the odd nonlinearity /(u) = |w|p ^u, for several values of p. 

(Hi) For p > 12, Dg is empty. 

(iv) For sufficiently large p, there exist /3 such that dec changes sign more than once as c varies from 
to c*. 
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